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Abstract 

In this paper we discuss the relation between the standard covariant quantum 
field theory and light-front field theory. We define covariant theory by its Feynman 
diagrams, whereas light-front field theory is defined in terms of light-cone time- 
ordered diagrams. A general algorithm is proposed that produces the latter from 
any Feynman diagram. The procedure is illustrated in several cases. Technical 
problems that occur in the light-front formulation and have no counterpart in the 
covariant formulation are identified and solved. 



1 Introduction 



Dirac's paper on forms of relativistic dynamics opened up a whole field of investigation: 
the study of different ways of quantizing and the relationship between different forms of 
dynamics. Three forms were identified, the instant form, that corresponds to ordinary 
time ordered theories, the point form, that will be of no concern to us here, and the front 
form, where the dynamical variables refer to physical conditions on a plane advancing with 
the velocity of light. (Such surfaces are called null planes or light fronts.) The latter form 
has the advantage that it requires only three dynamical operators, " Hamiltonians" , the 
other seven (kinematical) generators of the Poincare group containing no interaction. The 
other advantage noted by Dirac was that there is no square root in the Hamiltonians, thus 
avoiding the degeneracy of the free field solutions. (Particles and anti-particles cannot 
have the same kinematical momenta.) 

Dirac published his paper right in the middle of the period in the development of 
quantum electrodynamics, when old-fashioned perturbation theory was replaced by the 
covariant formalisms of Feynman, Schwinger and Tomonaga f^. Only much later were 
his ideas taken up againi. Quantum electrodynamics held center stage for a long time, 
not only for its unparalleled phenomenological success, but also because it functioned as a 
role model for many new theories, notably the gauge theories of the weak and the strong 
interactions. 

A new development occurred when the infinite-momentum frame, which had appeared 
in connection with current algebra (see e.g. De Alfaro et al. 0]), was proposed by Wein- 
berg P as a tool in the study of scalar theories, because it simplifies the vacuum structure 
in those theories. Not long after Weinberg's paper was published, it was suggested 
that the use of a new set of variables, viz, 

I ijC I ijC tjC "JO ]_ 2 / T \ 

would provide the advantages which are present in the use of the infinite-momentum 
frame. In such a description, plays the role of "time", i.e., the evolution parameter, 
and the connection to Dirac's front-form of dynamics seems immediate. (We will use the 
terminology light-cone time (Let.) for this variable.) However, the connection between 
the rest frame and the infinite-momentum frame involves a limiting procedure of Lorentz 
transformations. Therefore, the equivalence of descriptions in those frames cannot be 
derived using arguments based on Lorentz invariance alone. 

So the question concerning the relationship between different forms of dynamics re- 
mains difficult to answer. In particular, the connection between the manifestly covariant 
formulations and the front form is not yet fully clear. The main reason is that quantiza- 
tion using planes surfaces on which the initial conditions are specified-initial 
surfaces-is beset with difficulties, that occur already at the classical level in scalar theo- 
ries. It is a well established result from the theory of partial-differential equations |^ that 
the Cauchy problem with an initial surface that contains a light-like direction, is ill-posed. 

This disturbing fact might hinder the development of a Hamiltonian formulation of 
front-form field theory, if it could not be circumvented. A possible way out was shown 
by Chang and Ma and later by Kogut and Soper one may attack the problem 



^Leutwyler and Stern |^ extended Dirac's analysis to include two more forms of dynamics. A review 
on the present situation can be found in B. 



at the level of Feynman diagrams. If one follows that line, one must show that the 
usual Feynman rules can be reformulated in terms of the new variables, eq. (|T]), or their 
conjugate momenta 

In this paper we use this approach and show that it is successful for theories describing 
spinless particles. We show in sect, ^how to derive light-cone time-ordered (l.c.t.-ordered) 
diagrams from a given Feynman diagram by integrating over the I.e. energy p~ . The 
general algorithm is illustrated there by applying it to the box diaeram. On the way to 
the proof of equivalence we encounter questions of regular izationcl. For scalar theories 
they are not more difficult to answer than in the manifestly covariant formulation. 

The true difficulty lies in theories containing spinning particles. In the case of spin- 1/2 
particles one encounters the following expression for the free propagator |]10| 



p2 _ ^2 _|_ p2 _ ^2 _|_ 2p+ a P^ ~ + 

where pon is the on-shell value of the four momentum of the spin-1/2 particle with mass 
m, if its components p~^, p^ and are given. The component is computed from 

= 2p^p^ ~ p1 = 1^'^ (4) 

and so 

_ pI + 

Pan- • 

The occurrence of the non-propagating part i'-f^ /2p'^ makes the treatment of fermions 
in front-form field theories much more difficult, as it gives rise to integrals that are much 
more singular than the corresponding integrals in time-ordered or manifestly covariant 
formulations. In sect, ^the general case of spin-1/2 particles is discussed. As an illustra- 
tion a box diagram, describing two fermions exchanging scalar bosons, is reduced to a set 
of l.c.t.-ordered diagrams. 

The algorithm we propose demonstrates the equivalence of Feynman diagrams to sets 
of x"'"-ordered diagrams in the case of scalar and spin-1/2 fields. (In a way, this is the 
reverse of Wick's theorem for time-ordered perturbation theory. The fact that the Cauchy 
problem with a null-plane as initial surface is ill posed makes the Wick theorem in front- 
form dynamics a strictly formal result. Different interpretations have lead to different 
perturbative expansions [l^l-) The popular belief that massive fields do not have these 
problems is a misconception. The leading behavior of the fields near the light front is 
independent of the mass. 

A brief discussion on the extension of our treatment to diagrams with several loops is 
given in sect. ^ 

In the course of our investigation we encountered several technical difficulties. They 
are discussed in sect. ||, where solutions are given too. In particular we argue that there is 
no problem concerning zero modes, if the ^"-integrals are regularized properly, viz, using 
a regularization that preserves covariance. But it shows that Feynman diagrams give rise 
to terms in the perturbative expansion of the 5'-matrix that act on p'^ = 0-states. 



^The problem of renormalization is not discussed in this paper. 



The next section (sect. is concerned with the many mathematical details that were 
left out from the preceding sections, lest the main line of argument be blurred. 

We close with a discussion of our results and compare them to some of the literature 
on light-front field theory. 



2 Equivalence 

Before we proceed with the equivalence proof, we define what we mean by equivalence. 
By application of the Feynman rules as ordinarily understood, one obtains manifestly- 
covariant expression^ for terms in the perturbative expansion of S'-matrix elements, ex- 
pressed in terms of four-momenta, masses, spins, and dynamical ingredients: coupling con- 
stants. Wick's theorem can be understood as asserting that the S-matrix elements could 
be calculated as well in time-dependent perturbation theory, and as giving us an algo- 
rithm to combine the terms found in the latter case into manifestly-covariant expressions. 
Thus Wick's theorem establishes the equivalence of time-dependent ("old-fashioned") and 
covariant perturbation theory. 

In this paper we use the word equivalence in a similar way: each term in covariant 
perturbation theory-Feynman diagram-can be written as the sum of amplitudes that 
can be interpreted as terms in a Let. -ordered perturbation series. (In the interest of 
brevity, we will use the terminology l.c.t.-ordered diagrams.) In fact, those amplitudes 
are expressed in momentum-energy-space quantities, but it is a straightforward matter to 
translate them into space-time language, thus justifying our terminology. 

By taking Feynman diagrams as our point of departure, we avoid the problems of 
front-form quantization mentioned before. Besides, we also side-step the problem of 
identifying the independent degrees of freedom and the determination of commutation 
relations between them for a constrained system. 

The splitting of a Feynman diagram into l.c.t.-ordered ones results in amplitudes of 
the form 

Here, P~ plays the role of the "energy" variable, conjugate to the light-cone time (l.c.t.) 
x~^. Hq is again the energy, but now expressed in terms of the kinematical components 
and p± of the momenta of the particles in the intermediate state between two interactions. 
The objects V are the vertices that correspond to the local interactions. As we will show, 
expressions of this form arise naturally upon integration of a Feynman diagram over the 
the minus-component of the integration variable. In general, a number of l.c.t.-ordered 
diagrams are derived from a single Feynman diagram. This is directly analogous to old- 
fashioned perturbation theory, where n\ time-ordered diagrams sum up to one Feynman 
diagram with n vertices. However, there exists an important difference: in l.c.t.-ordered 
theory there are less diagrams owing to the linearity of the denominator of the single- 
particle propagator in the p~ variable 

1 1 , (7) 
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■^Except for non-covariant gauge terms in a non-covariant gauge like the light-cone gauge. 



Consequently, to every propagator there corresponds only one pole in p~ and its location 
in the complex p~-plane depends on the sign as well as the magnitude of p"*". This property, 
already alluded to by Dirac [Q, allows us to pose the condition that in any state p+ > 0. In 
the past the status of this so-called " spectrum condition" has remained somewhat unclear. 
We shall demonstrate that it follows directly from our splitting procedure and from a 
natural reinterpretation of amplitudes, quite similar to the reinterpretation of negative- 
energy states as states of positive energy of anti-particles in time-ordered perturbation 
theory. In fact, our reduction algorithm shows that the Let. -ordered diagrams have the 
property that the internal lines carry positive p"'"-momentum only. Therefore a better 
terminology might be spectrum property, but we stick to the term spectrum condition, 
because it is commonly used. In case one would like to formulate diagram rules in l.c.t.- 
ordered perturbation theory, one could use the spectrum condition as a limitation of all 
intermediate states to states where every particle has positive p"^-momentum. 

The spectrum condition is intimately related to causality. When the causal single- 
particle propagator eq. (|^) is Fourier transformed, one finds that the sign of p"*" determines 
whether one can extend the integral over the p~-axis to an integral along a closed contour 
in the complex p~-plane by adding a semi-circle at infinity for positive or negative '^p~: 
p"*" > corresponds to positive a:"'"-evolution. States with positive energy go forward in 
time and states with negative energy go backward in time. 

One can argue formally that the spectrum condition holds for all intermediate states. 
As a result of the completeness of the physical Hilbert-space each state is a superposition 
of free states and thus has a positive p"*" momentum. Any particle in a free state with 
positive p"*" has positive energy and goes forward in time. The conservation of kinematical 
momentum restricts the creation of particles in a Hamiltonian formulation, which 

is not the case in the equal-time formulation. We will show that this property indeed holds 
for the Let. -ordered perturbative expansion. 

This result seems rather obvious for spin-0 bosons, but, to our knowledge, its proof 
has never before been given. For spin-1/2 particles, there are complications due to the 
non-propagating part i'-y^ /2p^ of the fermion propagator. These render the equivalence 
proof in this case more difficult. 

There is an important point worth mentioning: the pole moves in the complex p~- 
plane as a function of p"*", and even crosses the real axis at infinity for p"*" = 0. This makes 
the propagator undefined as it stands. The crossing at infinity give rise to so-called 
"zero-modes" which will be dealt with later (see sect. |0| ). 

2.1 Examples 

As a pedagogical example we reduce the box diagram and the crossed-box diagram in 
0^-theoryi to the associated l.c.t.-ordered diagrams. This gives us the opportunity to 
show the working of the reduction algorithm in great detail. Later on we will give the 
general form of the algorithm. 

2.1.1 Box diagram 

The box diagram consists of four propagators: 



'^The type of theory is not essential for the arguments, the presence of a loop is. 
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Figure 1: Box and crossed-box diagrams. 
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with incoming momenta p and q and outgoing momenta / and p + q~l. The four momenta 
in the loop are: ki = k, k2 = k — I, k^, = k — p — q and k^ = k — p. We can rewrite 
the integral in terms of energy denominators and phase space factors. We define the 
quantities {Hi}f^^ as follows: 
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Then the k -integral can be written as follows: 

dr 

2^ {k- - Hi)(k- - H2)(k- - H3)(k- - H4)' 

The phase-space factor is given hy (f) — 16A;+(/c+ — l'^){k'^ — P"*" — Q'^){k'^ — P'^)- 

The positions of the poles in the complex k~ plane depend on the values of the external 
momenta and the value of k~^. To calculate the time ordered diagram wc must set the 
values of the external momenta p^, q^ and /+. For specific values of k^, the poles cross 
the axis and end up in the opposite half plane. When that happens, the integral changes 
discontinuously. The order in which the different poles cross the real axis depends on the 
values of the external momenta. In order to make our example definite, and without loss 
of generality, we assume p"*" > Z"*". Then we have five regions on the /c^-axis: 




(a) (b) 

Figure 2: The relative values of the longitudinal momenta, for the box diagram (a), and 
the crossed box diagram (b). All other box diagrams follow from these two through 
discrete symmetries. The absolute scale for the kf^s is set by the value of fc"*", the loop 
momentum. 



1. k+ < 0; '^Hi > 0, (1=1,2,3,4); 

2. < A;+ < /+; QH^ < 0, QHi > 0, (i=2,3,4); 

3. 1+ <k+ < p+; '^H2 < 0, '^Hi, S^ifs > 0; 

4. p+ <k+ <q+ + p+; QH-i > 0, QH, < 0, (i=l,2,4); 

5. p+ + q+ < k+; QHi < 0, (1=1,2,3,4). 

(see fig. ^) with respectively zero, one, two, three, and four poles below the real axis. Each 
region corresponds with a different process in light-cone time. This observation leads 
naturally to the definition of a skeleton graph. Each physical region in k~^ corresponds to 
one skeleton graph. It is a graph that is topologically equivalent to the original Feyman 
diagram, but has its internal lines graded + or — corresponding to the signs of the '^Hi, 
associated with the internal momenta fcf . 

In the first and the last case all the poles are at the same side of the real axis so the 
integral over k~ vanishes. Cases 2 and 4 are similar to each other. The integrals are 
calculated by closing the contour in the upper (lower) half-plane of complex /c~-values. 
The application of the residue theorem gives for case 2 



/^□ = — T^-TTZ T^—TZ r^, (ir 
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and for case 4 

/}4 ^ 1 I (12) 

° <j){H,~H,){H,-H,){H3-H2)' ^ ^ 

Case 3 is the most interesting one. Straightforward application of the residue theorem 
gives the result 

3 ^ (H, - H,){H, - Hs){H, + H,-H,- H^) 

° 4> {H,-H,){Hs-H^){H^-H2){H,~H,){H,-H2){H,-H2)' ^ ^ 




which can be spht into two parts 
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A point to be clarified is the meaning of the denominators {H^ — Hj). We choose 
the sign such that in these denominators > 0, corresponding to backward moving 
particles, while Hj has a negative imaginary part and refers therefore to forward moving 
particles. 

Four-momentum conservation gives = ki — p. We set k = ki and consider the case 
QHi < 0, QH4 > 0, which means that < A;+ < p+. Then, according to the residue 
theorem, we have a factor H4 — Hi in the denominator corresponding to this diagram. 
This factor can be written as follows 

(p± - k^Y + kl + m^ 
^'-^'=P - 2{p+-k^) W- ^^^^ 

The interpretation of the factor H4 — Hi is facilitated by cutting the diagram. We cut 
it first by a line cutting the legs of the loop with momenta ki and k^ and all incoming 
lines except p. For every internal line, we define an on-shell value of the corresponding 
minus-component as 



We define Hq as the sum of the on-shell minus-momenta on the lines cut. In our example 
we find for the internal lines 1, 4 and the external line q 
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The cutting line defines an intermediate state with total minus momentum P~ = 
p~ + q~ and total on-shell minus momentum Ho{l,4). The difference between these two 
is just H4 — Hi, eq. (|15|). Up till now, the direction of the internal four momenta is 
determined by the direction in which the loop is passed. If we reverse the direction of k^, 
the momentum with negative plus component, viz, k~^ — eqs. (|T^, ^), is replaced by 

— k^ , which is then correctly interpreted as the plus component of the momentum of 
the particle corresponding to line 4. 

If we consider a cut through ki, k^ and p, q, we will find k^ = k — p — q. The same 
algebra that led to the result eq. (\l7\} , will now give 
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It is clear that this procedure can be followed until the cut considered is cutting outgoing 
external lines only. There it stops. So, we conclude that we have the general result 
that any factor {Hi — Hj) is equal to the difference of the minus-component of the total 
momentum P~ = P + (1~ , and the on-shell minus-component of the momentum carried 
by the lines cut. Generally this holds for all combinations i and j such that QHi > and 
'^Hj < since the imaginary part is related to the sign of the on-shell momentum. 

The Let. -ordered box-diagrams can be interpreted as fourth order diagrams in l.c- 
perturbation theory, having the form 
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In order to systematize the reduction of Da to a sum of residues that correspond to 
Let. -ordered diagrams, it is appropriate to consider both the algebraic structure and the 
connection of residues with diagrams. 

First, we demonstrate the use of some concepts that will be of crucial importance for 
the proof of equivalence of the Feynman-diagram approach and Let. -ordered perturbation 
theory in the simple case of the box diagram. (The general proof is to be found in sect. 

The first object of interest is the Vandermonde determinant of order k: 
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The second one is Wn^miHi, . . . , if„|if„_|_i, . . . , -ff„+m) defined as 
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[21] 



By direct computation one verifies easily the following statements: 
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A(xi,...,x„,?/) nr=i(a;i -?/)' 

Straightforward application of our rules gives for the skeleton graphs the following corre- 
sponding amplitudes: 

^2 ^ i W,,,{H^,H,,H,\H,) ^ ^24) 
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In case 3 we have 
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needs to be rewritten such that energy denominators appear. The energy denominator 
(if4 — Hi)~^ should appear in all l.c.t.-ordered diagrams associated with D^. It can be 
extracted as follows: 
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and upon using eq. ( p2|) we recover the final expression (|j). The proof in sect, y demon- 
strates how this type of reduction can be carried out in the general case. 

Secondly, we describe the relation of this algebraic procedure with l.c.t.-ordered dia- 
grams. We begin with the Feynman-diagram and enumerate the possible configurations 
of poles in the complex A;^-plane. (Cases 1, . . . , 5.) A pictorial representation of those 
cases where the contour integral over k~ does not vanish (Cases 2, 3 and 4) is given by 
diagrams where the sign of ^H is indicated. (See fig. ^.) The box diagram is relatively 
easy to reduce to l.c.t.-ordered diagrams because in any associated skeleton graph there 
are at most two vertices that need to be ordered with respect to each other. There are 
only two internal lines which connect an incoming line with an outgoing line. Generally 
we call this kind of diagrams flat. The box diagram being flat, it does not show all the 
possible complications. Therefore, we also discuss the crossed box. 



(b) 




Figure 4: The crossed box. (a) Feynman diagram, (b) skeleton graphs , and (c) 
Let. -ordered diagrams. 

2.1.2 Crossed box diagram 

The difference between the ffat box and the crossed box is clearly visible in the sign 
patterns, for the imaginary parts of the denominators, one encounters when going around 

the loop. In the flat box one encounters the sign patterns H , + H and + + H — . 

In the crossed box, however, the sign patterns are H , H 1 — and + + H — . When 

there are two poles on either side of the real fc~axis, two sign changes occur, which can 
be seen as "bends" in the internal line, from backward to forward and vice versa. The 

skeleton graphs with signatures H and + + H — are treated in the same way as the 

corresponding fiat ones. The case H 1 — leads to four l.c.t.-orderings, which we explain 

now. There are two possible orderings of the two vertices with incoming external lines. 
Having chosen one ordering, we follow one of the internal lines until we reach a vertex 
with an outgoing external line. Either of the two vertices with outgoing lines can come 
first. This gives a total of four Let. -ordered diagrams. The corresponding diagrams are 

depicted in fig. ^. We choose again p+ > /+. The algebra for the H 1 — case, given by 

the interval — < < l^, gives the residue 
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that can be reduced to 



{H,-H,){H2-H,){H,-H3) ^ iH,-H,){H,-H3){H2-H,)- ^^^^ 

In order to expose the four l.c.t.-orderings, we write the two energy denominators next 
to the vertices with incoming hues as sums of two terms , e.g., 

1 1 
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Inserting these two denominators in the above expression gives four Let. -ordered di- 
agrams: 
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(31) 

The l.c.t-ordered diagrams and are obtained in a similar way as in the case of 
the box diagram. Now we can easily answer the question why the number of Let. -ordered 
diagrams is smaller than the number of ordinary time-ordered diagrams. In the latter 
case, a loop with 4 vertices leads to 4! = 24 diagrams. For the Let. -ordered diagrams 
the number is reduced, because of the smaller number of poles. We can also interpret it 
as the spectrum condition. The spectrum condition restricts the plus-component of any 
momentum on any line, internal as well as external, to non-negative values. Internal lines 
with negative plus-momentum correspond to poles with positive imaginary part and are 
interpreted as anti-particles. By reversing the direction of the momenta on these internal 
lines while maintaining four-momentum conservation, these lines can be again associated 
with particles. So, conservation of plus-momentum then provides the limitation of possible 
diagrams. The number of diagrams however, does depend on the external momenta; it is 
four in the case of the box diagram and six for the crossed box diagram. 



2.2 General case 

In a typical Feynman diagram on encounters several single-particle propagators. Following 
the same line of reasoning as used in the two examples given above, one will find that the 
corresponding poles in the loop variable k~ are located at different sides of the real fc~-axis, 
depending on the value of k~^. To illustrate this fully general property, we consider a one- 
loop diagram with vertices, N internal and external lines. The number of incoming 
lines is Ni; there are N2 = N — Ni outgoing lines. Suppose we call our integration variable 
k and identify it with k^ (see fig. |^). Four-momentum conservation takes the form 



(32) 




where the functions K-^{p'^, ...,p'^) are hnear. It is obvious that for arbitrary, but fixed 
external momenta all k^' < (fc/" > 0) for fc"*" — >• cxd (/c+ — >• —oo). So, we can divide the 
real A;"'"-axis into different regions, a semi- infinite region where all kf < 0, another where 
all kf > and — 1 finite regions where some of the kf are positive, the others being 
negative. 

We will use here again the concept of a skeleton graph, as we did in the cases of 
the box and crossed box diagrams. For any Feynman diagram with given values of the 
external momenta, and for every interval in the loop variable k~^, one draws a graph that 
is topologically equivalent to the original Feynman diagram and has its internal lines 
graded either + or — , corresponding to the signs of the imaginary part of the poles Hi. 
In the one-loop case we thus find for k~^ — > —oo the skeleton graph with all lines graded 
-|-; there are — 1 skeleton graphs with lines graded + as well as lines graded — and, 
finally, for fc"*" oo, a skeleton graph with all lines graded — . From our discussion of 
the causal single-particle propagator it becomes immediately clear that lines graded — 
(+) correspond to particles moving forward (backward) in x"'"-evolution. This justifies 
the terminology we adopt: if two vertices in a skeleton graph are connected by a line with 
internal momentum say k^' > 0, then the vertex from which the momentum /cf is flowing 
is said to be earlier than the vertex into which /cf is flowing. 

Apparently, each skeleton graph corresponds to a (partial) Let. -ordering of the vertices 
in a Feynman diagram. The graphs with all lines graded -|- or all graded — correspond to 
a cyclic l.c.t.-ordering of the vertices that contradicts logic. Fortunately, these graphs are 
associated with the situation that all poles in k~ are lying at one side of the real k~-a.xis, 
in which case the amplitude vanishes. In all other graphs there is at least one vertex with 
an outgoing internal line with positive plus-momentum and an incoming internal line with 
a negative plus-momentum, and at least one vertex where the situation is reversed. The 
former vertices are called early, the latter late vertices. A sign change in the skeleton 
graph corresponds with an early or a late vertex, the other vertices are called trivial, (see 
fig. P) If only one early and one late vertex are present in a given skeleton graph, the 
partial ordering is complete after the trivial vertices are ordered. This was the case for 




Figure 6: The sign change in correspond to early or late vertices. Both lines go in 
the same time direction. The sign of is opposite to the sign of kf. 



the flat box diagram. Otherwise, the different early vertices must be ordered with respect 
to each other and with respect to the late vertices. This additional ordering produces 
several Let. -ordered diagrams associated with a single skeleton graph. In this way, a 
single Feynman diagram gives rise to a number of consistently Let. -ordered diagrams. At 
this stage one reverses the directions of the four-momenta fcf on all those lines i where 
kf < 0. One sees immediately that as a result early vertices have only outgoing internal 
momenta, whereas late ones have only incoming internal momenta. 

We use the late vertices in a different way than the early vertices. Starting from an 
early vertex, all vertices on the lines going out from this vertex are ordered relatively to 
this vertex. However, vertices lying on different lines are not yet ordered relative to each 
other. For two lines starting at different early vertices and connected at a late vertex, we 
can fix the relative ordering of the intermediate vertices, since all vertices on both lines 
must occur before the late vertex. When all late vertices have been encountered, their 
relative ordering determines the complete ordering of the full diagram. 

The last point to be clarified is the form of the amplitudes corresponding to l.c.t.- 
ordered diagrams. We consider again the one-loop diagram and write the covariant am- 
plitude as follows: 

f d^k 1 

~ J {2ny[kl-ml + ze]---[k%-m% + ze] ^^^^ 

where for the sake of simplicity we have put all vertex functions equal to unity. Using 
eq. (p2D we can write for a typical factor in the denominator: 



kf - mi + ie = 2kt k^ - 



2kt \k~ + K. 



2kt J 
{k + Ki)l + m'i - ie" 



= 2kt{k--Hi) (34) 

The poles in k~ , Hj, are functions of the kinematical components of k^^ and the external 
momenta pj. The imaginary part of Hi, QHi, is determined by the sign of kf. Now 
suppose that for given external momenta k~^ is such that m pole positions are located in 
the upper half plane {QHi > 0) and n = — m in the lower half plane {QHj < 0). In 
order to simplify the discussion, we renumber the lines such that > for 1 < i < m. 



^Hj > for m+l<j<N = m + n. Consider the k -integral by itself: 

" / 2n2^ kt ■■■kt,[k--H,]--- [k- - Hn] ^^^^ 

After performing the integral by closing the contour in either Q'/c" > or Q'/c^ < one 
obtains a rational function of the His: 

Details on the functions Wm,n and A are given in sect. ^. As we argued before, the k~^ 
integration splits into intervals, (-oo, fc+(0)), (/c+(0), ■ ■ -, {k+{N -l),+oo), 

where the boundaries k^{i) are defined such that in the interval {k~^{r — 1), k~^{r)) there 
are, for finite r, precisely r iJj's with > 0. In each of the finite intervals one skeleton 
graph is present corresponding to one A;~-integral Dm,n- For either n or m different from 
1, Dm,n does not have the desired form eq. (^. The full Feynman diagram is recovered 
by summing Dm,n over m from 1 to iV — 1, integrating over k~^ over the appropiate finite 
interval, and over k±. A further reduction, corresponding to the transition from the 
partially ordered skeleton graphs to the completely ordered diagrams must be performed. 
The heuristics that help us to do so is provided by the space-time concepts. Take any 
early vertex and identify it with an event at some Let. = tq. Suppose Hi and Hm+i 
are the poles corresponding to the outgoing and incoming internal lines resp. at this 
vertex. The intermediate state with momenta ki and fcm+i corresponds to the l.c.-energy 
denominator Hi — Hm+i- We can convince ourselves that this is correct if the reversal of 
four momenta on lines with fc^ < is effected. We find 

Hi - H^^i = P - - ^"^+1^ + . (37) 



2k I 2k, 



m4 



The identification of becomes clear when one uses the surface r to cut the internal 
and external lines. We can close the surface by cutting all the external lines attached to 
the piece of the Feynman diagram out of which the p"*" momentum flows. The momentum 
flow through this first cut plane Tq equals the flow through the second cut plane r since 
energy-momentum is conserved locally in a Feynman diagram. The flow of P~ through the 
internal lines equals the initial flow into the diagram minus the flow through the external 
lines cut at r. The loop momenta do not contribute because they go into the r cut plane as 
well as out of the cut plane, so there is no net contribution of these momenta. (We stress 
here that this interpretation is correct only after reversal of the four momenta on the lines 
with QHj > 0.) For brevity we call denominators of the form H^ — Hj = P~ — HQ{i,j) 
energy denominators. These cut planes can be interpreted as equal-time surfaces. 

Now the strategy is clear. For every skeleton graph one uses the surface x+ = r to 
cut lines that give rise to energy denominators. That this is possible is the content of our 
proof of equivalence. Indeed, as we demonstrate in sect. 13, we have 
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Figure 7: Illustration of the cut-procedure in a general n-leg diagram. For each cut 
area there is four-momentum conservation in a Feynman diagram. We can use this to 
make general statements for an intermediate state. The minus-momentum transferred 
across the intermediate state is equal to the incoming minus-momentum. The sum of 
energy denominators contains the sum of all external minus-momenta minus the sum of 
all on-shell values of the particle lines. Loop momenta drop out because these go in and 
out of the cut area. 
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Figure 8: The product of two time orderings is the sum of the relative orderings. 



That this identity leads to a recurrence follows from the fact that the two terms in the 
r.h.s. of eq.(P5|) have the same form as the original one. The reduction stops if either n 
or m is reduced to 1. 

There remains one loose end that we tie up now. If several early vertices occur, we 
have to consider all l.c.t.-orderings of them. (This happens if n > 2.) Moreover, we 
may have to consider l.c.t.-orderings such that some late vertices occur before some, but 
not all, early vertices. In all those cases the number of internal lines cut by an x"*" = r 
surface is greater than two, but always even: for every line going into this surface there 
is a corresponding outgoing line. We have seen that a pair of lines, one incoming, the 
other outgoing, that connect in an early (or late, for that matter) vertex, gives rise to 
an energy denominator, say Hi — Hj. When two pairs of such lines occur, there will be 
two energy denominators, which we call simultaneous parts. A simple example illustrates 
this. Let the two early vertices be a and and and pj^ the momenta on the two 
corresponding incoming external lines. The reduction algorithm gives two factors, one 
corresponding to the vertex a, of the form l/{Hi^ — Hj^), the other being l/{Hi^ — Hj^). 
As before, we can rewrite such factors, upon reversing the backward flowing momenta, in 
the form l/(P-(/5) - Ho{(3)) and 1/{P (a) - Ho{a)), where p-(7),7 = a, (3 is the total 
net external minus momentum flowing into vertex 7. A simple algebraic identity 

(39) 



P-{a)-Ho{a) P'{P)-Ho{P) 



1 /I 1 



P~{a) + P-{P) - Ho{a) - HoiP) \P-{a) - H^{a) P~{P) - Ho{P) J 

combines the two factors in the correct way. The first factor can be rewritten as l/(P^(aU 
P) — Hola U P)), which we recognize as an energy denominator for the intermediate state 
with the four internal lines ia, ja, i/s and j/3. The combination l/(P~(a) — Ho{a)){P~ {aU 
P) — H^i^aU P)) corresponds to the l.c.t.-ordered diagram where the vertex a comes before 
vertex /3, the other part of the r.h.s. of eq.(^) corresponding of course to vertex P 
preceding vertex a (see fig. |). Clearly, the splitting formula works also recursively, so it 
applies to any number of pairs of internal lines. 

With this observation we end our general discussion of the reduction algorithm. We 
stress that the l.c.t.-ordered language used here has heuristic value only, but does not 
replace a strict proof. The algebraic details are provided in section ^. 



3 Spin-1/2 particles 



In the previous sections we dealt with scalar particles only, thus avoiding complications 
due to summations over spin degrees of freedom in intermediate states. Here, we discuss 
these complications for spin-1/2 particles. The reduction algorithm in this case is partly 
identical to the algorithm for scalar particles. However, we now have to include in our 
treatment not only the energy denominators, but also the numerators. Consider the 
Feynman propagator for a single particle. The spin sum + m depends on so we 
have to account for that when we define the skeleton graphs corresponding to a Feynman 
diagram. It has been argued before that one can split the Feynman propagator into 
two pieces, one that is independent of p~, the instantaneous part, and another piece, the 
propagating part, where p~ occurs in the denominator only: 

p2 _ ^2 _|_ p2 _ _|_ 2p+ ^ p^ — + ie 2p+ ' 

with the obvious definition p^„_,^^,; = {p'on-sheii^V' ^Vi).^^^"^^ V'on- shell = {pl+m^)/{2p+). 
The spin sum J2a'^^°'Kp) '^u^^'^p) runs over a complete basis in spin space, viz, particles 
and anti-particles. 

In order to illustrate the differences between the purely scalar case and the situation 
where spin-1/2 particles occur, we discuss in the next subsection the flat box diagram 
with two bosons and two fermions. 



3.1 Example: fermion box diagram 

Consider the case of two spin-1/2 fermions exchanging spinless bosons. A Feynman 
diagram for the fourth order in perturbation theory is shown in fig. (a). The momenta 
are defined similar to the scalar case, fig. Before the associated skeleton graph can 
be drawn, one must split the fermion propagators for the intermediate states into the 
two parts: instantaneous and propagating. This results in four different diagrams shown 
in fig. ^ (b). In these four, the ^'-dependence of the propagators is present in the de- 
nominators only. Therefore, one can apply the methods described in the previous section 
immediately, since the numerator does not depend on the integration variable. Doing 
so, one obtains the eight skeleton graphs drawn in panel (c). These graphs form the 
basis of the splitting into l.c.t-ordered diagrams. After this has been achieved one can 
combine certain diagrams into a single diagram by adding the propagating part to the 
instantaneous part of the fermion propagator. This regrouping of diagrams is represented 
graphically in the last panel of this figure. The formulae associated with the four final 
diagrams are 
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Figure 9: (a) The Feynman diagram, (b) the corresponding diagrams with on-shell spinor 
projections or instantaneous parts, (c) the skeleton graphs, (only the tilted box leads 
to two l.c.t.-ordered diagrams) (d) the summed Let-ordered diagrams, which yields the 
adjusted blink propagators. 



(27r)30 (P- - iJo(l, 4))(P- - i/o(l, 3))(P- - Ho{2, 3)) ' 
(27r)30 (P- - Ho{l, 4))(P- - i7o(2, 4))(P- - Ho{2, 3)) ' 



^ /■ dfc+d^A;^ 7°A37'' (S) YMl'^ . . 
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The objects A and Q are defined as follows: 

Ai = hon + rrii, (42) 

and 

a = ^i + m, = Ai + 7+(A;r - A;-). (43) 

The on-shell values of the minus components have been defined before, see eq. (^). The 
energy denominator P~ — Hq{1, 2) is of course 

P- - Ho{l, 2) = p- + q- - k^^^ - k^^^^ - 

2k^ 2A/2~ 

and the other ones are defined similarly. The phase-space element (f) has also been given 
before, below eq. (ITO). 



The fiat fermion box shows clearly the peculiar complications caused by spin. Because 
the numerators and the denominators of the fermion propagators depend both linearly on 
the integration variable k~ , one has to perform a Laurent expansion in order to identify 
the pole terms. This leads to a number of "intermediate amplitudes", equal to 2^, where 
F is the number of internal fermion lines. These amplitudes give rise to skeleton graphs 
that can be reduced to l.c.t.-ordered diagrams in the, by now, familiar way. The l.c.t.- 
ordered diagrams in which an element --y^ /2kf occurs, are special, because associated 
to every one of them there occurs a diagram where the element '-j^ /2kf is replaced by 
h.i/{2kf{P^ — Hq)). This happens only in those cases where the on-shell value Hi = 
{kf + mi)/2kf occurs in a single energy denominator. These are the states that begin 
and end with the creation and the annihilation of the same particle. We call an internal 
line with this property a hlink. In the case considered above: D*-^^ contains in only 
one factor in the denominator, the same holds for D^^^ and i/i, whereas D^^-* contains the 
blinks Hi and H^. If a blink occurs, one can recombine, after the Let. -ordering has been 
performed, the propagating part and the instantaneous part into a complete propagator. 
This is done in the diagrams P)(i), P)(2) and P)(3) gq_ (|2|), and illustrated in fig. |(d), 
where the thick lines beginning and terminating in dots symbolize complete propagators. 

3.2 Including the instantaneous terms 

By now, it is relatively easy to formulate the general reduction algorithm. It has four 
steps. 



(i) For a given Feynman diagram, perform the Laurent expansion of the fermion prop- 
agator, i.e., split the propagator into a propagating part and an instantaneous part; 



Figure 10: The state begins and ends with the creation and annihilation of the same 
particle. The singularity in this diagram cancels against the same singularity in the 
instantaneous diagram. 



(ii) Determine the skeleton graphs for all diagrams obtained in step (i); 

(iii) Perform the reduction of all skeleton graphs in exactly the same way as it was done 
in the scalar case; 

(iv) Identify the blinks and sum the diagrams corresponding to the same blink in order 
to obtain amplitudes with complete spin sums. 

In order to understand why we recommend step (iv) we consider the general case. Let 
/cf be the four momentum of a blink. The two corresponding diagrams, partly shown in 
fig. |10|, contain the factors! 

1 A. 



2ktp--{---Hr--y 

Ai = -f^^Hi + -i'kf + 7^ ■ A;^ + m^, 

Qi = i-'/2kt. (45) 

We see that both Gi and Qi become singular at fc,^ = 0. These two singularities appear to 
cancel. We can see this most clearly if we realize that the denominator P~ — {■ ■ ■ Hi ■ ■ ■) 
can be rewritten as k^ — Hi. If kf — >■ 0, then Hi oo, so we see that in this limit Gi 
behaves as follows: 

Therefore, the sum of the two contributions is finite for kf ^0. 

_ 1 l~^Hi + -^-kf + 7_L ■ A;± + mi 7+ 
'^1+ Qi — 777+ T-z 77 r 



2kt k- - Hi 2kt 

j+k' + 'j~kt + 'y±- k^ + rrii 
2kt{k--H{) 

'j+k' + 7_L • fc± + 



(47) 



^We define 7_l = (0,71,72,0), k± = (0,fci,fc2,0) and ■ kj_ = -(71^1 +72^:2)- 



b 

Figure 11: The two- loop diagram. 



This expression appears after the integration of the energy variable. Then k~ is a function 
of the external variables only and represents the p~-flow through the intermediate state 
under consideration. 

3.3 General case 

We have seen that singularities in l.c.t-ordered diagrams cancel similar singularities in 
instantaneous terms. The instantaneous terms might contain other divergences that are 
cancelled by lower-order terms with additional instantaneous terms. The question re- 
mains whether this procedure ends, or whether we are left with terms which contain only 
instantaneous singularities that do not cancel each other. In the section on divergent 
contour integration we show that the proper treatment of the shift of poles to infinity 
removes all singularities from each residue. So after the recombination of terms we won't 
have a residual term in the form of instantaneous parts. 

Although we are not concerned with gauge theories explicitly, we note that most of these 
terms drop in a gauge theory with the (naive) light-cone gauge, and in theories with scalar 
and pseudo-scalar couphng, due to 7+7+ = 7+7^7+ = = 7+7*7"*". 

4 Multi-loop diagrams 

The extension of the reduction algorithm from Feynman diagrams with one loop to Feyn- 
man diagrams with several loops is not difficult, but there are some points that need to 
be clarified. The loop integrations can be done one after another since the structure of 
a l.c.t.-ordered diagram is not essentially different from a Feynman diagram. We will 
illustrate the procedure with a simple example in the section below. 



4.1 Two- loop diagram 

We consider the scalar diagram with two loops, depicted in fig. |ll|. The corresponding 
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Figure 12: The imaginary signs of {Hi, H2, H^, H4, H^} for the different sectors in k'^ 
space. Only the inner regions I and H correspond to integrals and skeleton graphs. 



integral is : 
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dg d/c 
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Where the phase factor = 2^k^{k^ — p^){k^ + q^)q^{q^ and the poles are given 

by: 
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(49) 



2q^ 



There are twelve sectors in k'^ g^-space corresponding to twelve skeleton graphs. These 



sectors are depicted in fig. where also the signatures of the skeleton graphs are shown. 
The amplitude D vanishes if either the integral over k~ or the one over q~ vanishes. 
The former happens if Q^i^i, 55-^2 and '^H^ have the same sign, the latter if this happens 



for Sif3,9=if4 and ^Hr^. We read off from fig. 12 that there are two sectors remaining, 



denoted as I and H. In sector I we have Hi, H^, H4 < and H2, H^ > 0. The reduction 
algorithm gives the l.c.t.-ordered diagram I of fig. |13|, with Let. ordering a < b < c < d. 
In sector H we have Hi^H^ < and H2, H^, H^ > 0. The corresponding Let. ordering 
is a < c < b < d. The only difference between the two diagrams is the sign of QH3, 
that is linked to the two different Let. orderings of the vertices b and c. For the sake of 
completeness we give the algebraic expressions for the two Let. ordered diagrams. Upon 





Figure 13: The two l.c.t-ordered diagrams that follow from the two-loop diagram. 

integration over k~ we obtain: 

i f dq~ 1 



271 J <P (H, - H,){q- + H, - H,){q- - H,){q- - H, 



p, ^ /■ 1 

" = 2^I^{H,- H,){q- + H,- H,){q- - H,){q- - H,)' 

The q~ integration is straightforward and gives the result: 

1 1 

Di = 



_ 1 1 
" " {H, - H,){-H, -H, + H,){H, - H,) ^^^^ 

After reversing the directions of the lines corresponding to negative kf , viz, k2 = k — p 
and fcs = g in Dj and k2, k^ = k + q and A;4 = p + g in Dji we obtain the Let. -ordered 
diagrams depicted in fig. 0. As before, every factor in the denominators of eq. (^) can 



again be written in the form P~ — HQ{i,j), where HQ{i,j) is the sum of the energies on 
the lines i and j. 

We see again that the integrations over k~^ and q~^ are limited to finite regions. After 
reversing the lines with negative QHi, one sees that the diagrams obtained have the 
spectrum property. 



4.2 General multi-loop diagrams 

In an arbitrary Feynman diagram with L loops, one must first identify the independent 
integration variables, say gf,...,g^. Then one can characterize the different types of 
pole positions in an L-dimensional space with coordinates {qt,---,QL)- The different 
signatures C^Hi, . . . , '^H^) divide this space into a number of sectors, each sector being 
associated with its particular skeleton graph. The sectors where the pole positions in 
all variables q~ are distributed over both half planes, Sg," > and '^q^ < resp., are 
necessarily finite. This is so, because for any loop i, all poles Hj^ occurring in this loop 
will have '^Hk < {QHk > 0) if the integration variable q^ goes to infinity (—infinity). 
Therefore, the sectors in [qf, . . . , g^)-space which are semi-infinite in either of the qf do 
not contribute to at least one of the integrals over the g,~ variables. 

So, in general we will have a finite number of skeleton graphs that each gives rise to 
a finite number of Let. -ordered diagrams. Each and every one of them has the spectrum 
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Figure 14: Two cases where the energy integral is ambiguous: (a) If the momentum 
along the loop is constant, (b) If there is one pole left in the the ^"-integration; the pole 
of the boson accompanied by the instantaneous part of the fermion propagator. 



property. In the case of spin- 1/2 particles, one can duplicate this algorithm, provided 
the full Feynman diagram, containing F fermion lines, is first split into 2^ intermediate 
diagrams according to the division of the spin-1/2 propagator into instantaneous and 
propagating parts. Then the reduction algorithm is applied to each of the intermediate 
diagrams, giving rise to the appropriate skeleton graphs and finally to the Let. -ordered 
diagrams, as was demonstrated in the one-loop case in the previous subsection. Of course, 
in the multi-loop case blinks may occur as well as in the one-loop case. They are treated 
in exactly the same way as before. Thus we see that the multi-loop Feyman diagrams, 
although algebraically more involved than the one-loop cases, can be reduced to l.c.t.- 
ordered diagrams using precisely the same algorithm as was used for one-loop Feyman 
diagrams. 

A final remark concerning the ie-prescription is in order here. It is used to define the 
deformed integration contours in all variables q~ simultaneously. After the residue theo- 
rem is applied to perform the contour integrals, the real parts of the poles are substituted 
in formulae like eqs. ( |50| , pl| ) (e = 0). If one would substitute complex poles in eq. |50| , 
ambiguities might arise in the values of the imaginary parts of the poles in q~. 



5 Technical difficulties 

In the previous sections we dealt with the equivalence between Feynman diagrams and 
l.c.t-ordered diagrams when the integration over k~ is well-defined. There are two types 
of special cases were the k~ integration is not well-defined. We can best illustrate these 
with simple examples. Consider a scalar loop like in 0^ theory. If the external line has 
positive p'^ then the integration domain in k~^ is the interval < k^ < p'^ (see fig. |l4 
(a)). One may wonder what will happen if p'^ = 0, because in that case the measure 
of the integration interval is zero. The poles in the two propagators cross the real axis 
in the /c~-plane for the same value of k^: k~^ = 0. If this diagram is finite, there must 
occur a delta-like contribution in k^. In cases where < 0, one can approach p"*" = 
by performing a Lorentz-transformation (that, however, does not belong to the stability 
group of the null plane) and take the limit. Such a Lorentz-transformation is always 
possible for a space-like external momentum, and there are situations where the momenta 
on three space-like external lines can be transformed to have = simultaneously. 
In other cases, like in (generalized) tadpoles, the measure of the integration domain is 



rigorously zero. We will consider a general approach which holds in all cases and give 
the same answer as a covariant calculation in the limit = 0. Tadpoles have a close 
relation with the ordering of operators in the Hamiltonian, therefore we see that the 6{k~^) 
contributions have a relation with the ordering. 

The other case where the k~ integration is ill-defined, occurs if at most one pole is 
present in the k~ integration. This happens for all Feynman diagrams with at most one 
boson propagator and at least one fermion propagator in the loop. Then the intermediate 
diagrams with an instantaneous part of the fermion propagator combined with a boson 
propagator needs regularization. Other examples are diagrams with at least two instan- 
taneous terms, which also lead to divergent integrals. The first order correction to the 
fermion self energy in a theory of scalar bosons and fermions with Yukawa coupling, is a 
simple example (see fig. (b)). The fermion propagator has an instantaneous part such 
that the only ^"-dependence of the integrand resides in the boson propagator. This inte- 
gral is not defined, so we need a way to deal with this type of integrals in a consistent way. 
We are primarily interested in a treatment which does not interfere with the algebraic 
rules, so the regularization must be a linear operation. In addition we require it to be 
homogeneous in the integration variables. In both cases, one where tadpoles are present 
and the other where instantaneous parts give rise to infinities, we are lead by covariance 
in our choice of regularization. Other arguments do not restrict the regularization to a 
unique method, while covariance does. 

5.1 "Zero modes" from energy integration. 

One of the integrals which show the presence of zero modes in a time-ordered formulation 
has been discussed already by Yan : 

f ^ 1 
lim / dp = 

k+^oJ (2(p+ — k+)p^ — ra^ + ie){2p+p^ — ra^ + it) 

1 — 1 covariant iTT 5 (^p ) / r r\\ 

/o + rr~^ = — ~ ■ (^2) 

[Zp+p — + zej"^ 

For p+ 7^ there is one double pole either above or below the real axis so it was concluded 
that, since p+ = is an unphysical value (no free state can acquire p"*" = 0), the integral 
should vanish. A careful analysis shows that eq. (|52D is an ambiguous expression, so one 
can get any value (including the covariant answer) and one needs to choose a regularization 
to get a well-defined integral. (The proper covariant value was obtained by Yan by taking 
the limits p"*" | and p"*" | in a special way.) 

For p+ = the p~ integral diverges, so (p+ — > 0,]?^ ^ oo) is the ambiguous point. 
If p"*" moves along the real axis and crosses p"*" = 0, the poles move through infinity and 
end up on the other side of the real axis. To deal with all singularities of this type at the 
same time, we introduce the variable u = l/p~ and study a general case: 



Dn = J dp- 



1 



i2ptp- - Ei + %e){2ptp- - + 2e) ■ ■ ■ (2p+p- - E^ + le) 

(2p\ - [Ei - te)u) ■ ■ ■ (2pt - m - ^e)u) ■ ^ ' 



The integrand goes to zero like for u — >■ ±cx), therefore the integral is well-defined, 
unless the integrand has a singularity at u = 0. So the only divergence can occur if u — >■ 
which gives a finite contribution only if all p^-momenta vanish at the same time. The 
poles in the variable p~ that moved to infinity, now correspond to poles in the variable u 
that cross the real axis at m = when either of the variables pf is zero. If all pf happen 
to be equal, the integrand is singular at m = 0. (The first example, eq. (^2|), is the special 
case with n = 2;H^ = H2 = m?.) This gives a finite contribution to the integral of D 
over on support = 0, thus D contains a delta function in 

The u coordinate regularization replaces all other arguments we might have to deal 
with this "zero- mode" problem. The choice of regularization determines the integral 
uniquely. Instead of treating the general expression, eq. (|53|), we regularize the case of a 
single pole in the integrand and use an algebraic relation to obtain the general expression. 
Consider the integral 

Di= j du / . (54) 

This expression is ambiguous for two reasons: it has a pole at m = and a double pole 
occurring for p"*" = Am = 0. The first ambiguity we remove by adding a small imaginary 
part to one factor u coming from the Jacobian. In order to obtain a covariant expression 
we must do this symmetrically: 

^ ^' ^ +^)- (55) 



u 2 \u + i5 u — iS 

We split the integral into two pieces; one just above the real axis and the other just 
below it. We do not give the singularities some strict nature (like principal value), which 
would lead to the square of the principal value for Generally we treat the energy 
as a complex variable (since each pole corresponds to a particle), and the kinematical 
variables are treated geometrically. The choice eq. (^) separates p^ > from p"*" < for 
all positive values of e and 6. Thus we find for the regularized integral: 



1 



2\u + i6 u - i5 ) {2p+ - {H^ - ie)u) 
ni9{p^) TTi9{—p^) 



+ 



2p+ + i{H^ -ie)6 2p+ - i{H^ - ie)6 
2m9{p+) 



2p+ + i{H^ -te)6' 



(56) 



For < we reversed its sign to obtain the last line. Integrating D^^^ over p"*" from 
to a cutoff A and taking the limit e — >• 0, gives 7r{\n{H-^) +\nS + 7ri/2 — ln2A). We 
shall see that the constant part ln(5 + 27r/2 — ln2A drops if we have two or more energy 
denominators. Using the algebraic relation 

1 " 1 

fn^+^- rr± , • Nrr 7T71 (57) 



the regularized integral becomes: 



^dp- 



1 ^ m5{p+)\n{Hi 



T.Tr^9wr^- (5^ 



i2p+p- - Ht + ze) • • ■ {2p+p- - + le) ti Tl.AHt ' H,^ 



The function S{p'^) appears here, because the integral is strictly zero for p"*" 7^ 0, although 
integration over p"*" gives a finite result. The result eq. (0) can also be obtained as 
a limiting case of eq. (|53D in the simultaneous limits pf —>■ One can check that 

eq. (|58| ) is the same as the covariant result, using a Wick rotation such that 2p'^p~ —p\^ 
— The limit for Hi Hj is well-defined. The constant term has dropped since 

k=i Uj^kiH,, - Hj ) 

which follows from the fundamental theorem of algebra. 

We emphasize here that these zero modes appear in loops where the p"*" momentum 
is constant along the loop. Zero-modes can be interpreted as an infinite number of states 
(around p^ = 0) which are infinitely off-shell = 00), and thus have zero probability 
for propagation over a finite distance. The combination of both gives a finite contribution. 
This is reminiscent of the ultra-violet divergences, where the large number of high-energy 
states give an infinite contribution. 

Although zero-modes are needed to obtain the covariant answer, they remain slightly 
artificial, which can be seen from the configurations were they occur. It seems that nature 
is telling us that the high density of states for high energy causes trouble: divergent inte- 
grals appear which have to be regularized, and zero-modes. Both result from singularities 
on the light-cone. 



(59) 



5.2 Divergences in the fermion loop 

In the section on diagrams containing fermions we stated that they could be reduced to 
Let. -ordered diagrams, provided no additional singularities would occur. In this section 
we deal with these singularities. We state that the regularization proposed here removes 
all of them. There remains one point to clarify, i.e., whether the method of regularization 
does indeed produce a covariant result. The latter point, however, will not be discussed 
in this paper. 

We now have the tools to deal with the singularities in the fermion loop. Earlier (sect. 
^ we saw that a blink combines two singular terms in such a way that we get a non- 
singular expression. However, in general the low-order terms, with several instantaneous 
contributions, are singular by themselves. We will show that the contribution from the 
contour at infinity leads to these singularities. After subtraction of the latter the singu- 
larities are gone and a proper recombination of terms will remove apparent singularities. 

For Feynman diagrams with at most one boson propagator in a loop there are singular 
parts in the contour integration. Even if the diagram would be convergent in the ordinary 
sense, i.e., in the covariant or instant-form, it will still be divergent in ^"-integration. 
The singular behavior of the fermion propagator on the light front is to blame. 



Also in the case of bosons an ambiguity occurred (see sect. |5.1|) . The result did 
depend on whether the contour was closed in the lower half plane or in the upper half 
plane. This problem could be resolved by choosing a particular combination of contours. 
We considered two contours, one consisting of the real axis and a semicircle in the upper 
half plane, the other one has the semicircle in the lower half plane. The integral over the 
real k~-axis was replaced by the average of the integrals over these two contours. This 
regularization turned out to have a number of desirable properties. 



Figure 15: The difference between the upper-plane contour and the lower-plane contour 
is the contour at infinity Too- 



In the case of fermion loops the problem is more complex. A straightforward appli- 
cation of the residue theorem gives results that depends strongly on the choice of the 
contour. The integral is also more divergent than one would expect from naive power 
counting, and has non-covariant singularities. At the origin of these problems lies the 
contribution of the (semi-) circle at infinity to the contour integral. We have to subtract 
this contribution. In a sense we propose a regularization of the contour integral. 

For an integral that converges on the real axis, it does not matter whether we close 
the contour in the upper or the lower half-plane. This means that the sum of all residues 
is zero: 

Resi = 0. (60) 

i 

If the integration is divergent we have contributions at infinity, which add differently to 
the two contours (see fig. ^). Therefore the two contours give different results. We 
denote this difference by Foo- Then we have 

^ReSi = Foo. (61) 

i 

To solve this problem we "regularize" the residues by subtracting fractions of Fqo from 
them: 

Resr^ = ^(Res, - a,T^) = 0, = 1. (62) 

i i i 

We then find the desired result that closing the contour in the upper half plane gives the 
same result as closing it in the lower half plane. We will see that the a^'s can be chosen 
such that the singularities cancel. 

The /c~-integral corresponding to a general Feynman diagram with a fermion loop can 
be written as: 



J 2-kt---k+ {k- -Hr)---{k' -H^) ■ ^^^^ 

We wrote only the structure that depends on k~ and k^ explicitly. We can consider one 
residue at a time, since the singularities of one residue are not canceled by another residue. 
The residue of the pole {k~ — Hj)~^, and the contribution of the contour at infinity are: 



I (7+^, + A)---(7+g,+/ 3n) 

2'^kt---k+ iH,-Hi)---{H,-H^ 



Res, = ,^^'\rr rri )rr rr (64) 



Too = ^ . ^+ ((E g.) (7-^ • • • 7+) + (A • • • 7+) + • • • + (7-^ • • • Pn)j (65) 



The latter can best be calculated by changing the integration variable to m = {k~)~^ and 
integrating over a circle around the origin with a small radius, |m| = e. The residue ReSj 
has only singularities if fc^ — 0, since then Hj oo. In the limit kj'—>-0 only a few 
terms survive. The surviving structure is precisely the contribution from the contour at 
infinity, which is independent of j 

lim Resj = Too • (66) 

J 

If we decompose the contribution from the contour at infinity now in such a way that 
each term contains only one singular point [kj' 0), we can subtract these terms from 
the residues with the same singularities, and the resulting regularized residues are finite. 
If we define the quantities and as follows: kj~ = k~^ — and Hj- = 2kj~Hj, then 
we find that the following regularized residues are finite for kj^O 

Res^^ = ReSj — ( ^„ —Hj + 



2^{k+-ql)U{qt-qty 



'^l ^3 \ f -\- -\-\ 

\ 2-+Hfc+ - - 4) n(g; - qt) ^ nqt - qt) mt - ^ '"^ 

1 




Since each residue is now regular we know that there is a combination of Let. -ordered 
diagrams where each term is convergent. The blink procedure tells us how to do this 
step-by-step. The argument above tells us that there are no singularities left. 

We killed two birds with one stone: we resolved the ambiguity in the contour inte- 
gration and removed the singularities. These two problems are intimately related since 
the singularity comes from the pinching of the contour at infinity between the "pole at 
infinity" and an ordinary pole. 

Another advantage of this procedure is that only the physical sectors are non-zero, a 
property that would be destroyed by ordinary contour integration. It allows us to keep 
a simpler view of causality and unitarity, where each line is associated with a particle 
moving forward in I.e. -time. In non-physical sectors all the poles are on the same side of 
the axis. Then formula (^) tells us that the result is zero. We do not know whether this 
regularization leads to the same amplitude as the covariant calculations. 



5.2.1 Example; vacuum polarization 

We will illustrate the procedure from the previous section with an example. The simplest 
diagram is the vacuum polarization diagram; a closed fermion loop with two fermion 
propagators. 

J.,. ^ f dfc- Trlril^k- + f3,)Y{l^k- + /j^)] 

J 2m Akfktik- - Hi){k- - H2) ' ^ ^ 

We deal with the physical sector so we can assume, without loss of generality '^Hi > 
and 53i^2 < 0. The result depends on the way the contour is closed. We will close the 
contour in the upper half plane. In terms of Let. -ordered diagrams we have the ordinary 
diagram, with two propagating fermions, and the diagram with the instantaneous part 



associated with the pole H2 in the lower half plane. Together they give the residue of the 
pole Hi, as expected. 



If we would have chosen to close the contour in the lower half plane the result would be 
different. The difference is the result of a finite contribution of the semicircles. We still 
have to subtract the fraction of the contour at infinity which is given by 

°° Aktikt - kt) \ ' ^ k^ - kt + kt - kt) ^ Aktikt - kt) ^ Aktikt - ktr 

(70) 

The l.c.t.-ordered diagram obtained has no singularities, and is independent of the direc- 
tion in which the contour is closed, since it is symmetric in Hi and H2: 



Tr[7n( ^-frg^- ) 7^ + /5i)7^( (^^f5|^) 7+ + 
Ak+kUHi-H^) 



or - aiT^ = D^:^ = ^ ' ^ ' ' • (71) 



This result could not be obtained if we would have taken a combination of an upper-plane 
semicircle and a lower-plane semicircle, this would give a singular, and thus an ambiguous, 
result. The contribution of the contour at infinity should be decomposed in a unique way 
to obtain a regular expression. 

We will not calculate the amplitude here since the diagram is divergent and the comparison 
is with other results is spoiled by renormalization. After this regularization the fc^- 
integration is automatically finite, since the domain is finite and the integrand is regular. 



5.3 Dynamical spin 

In a Hamiltonian theory we separate kinematical variables from dynamical ones. Starting 
from one equal-l.c.t surface, we evolve in the time direction to the next equal-l.c.t. surface. 
In a covariant formulation, different interaction points could be on the same equal-l.c.t 
surface. Constraints should, in a Hamiltonian formulation, account for these parts. A 
straightforward interpretation does not exist, the constraints "evolve" with the order in 
perturbation theory, and so does the physical Hilbert-space (constraints are working on 
the Hilbert space). We will illustrate this point with a Hilbert-space interpretation of the 
instantaneous interaction '-j^ /2p^ . 

The (p^)~^ singularity is ambiguous as it stands. We need a way to look at this 
part such that we avoid additional divergences which seem to appear. Note that the 
spinor-matrix elements of 7"*" are the same as those of . Therefore we might suspect 
that the occurrence of 7"'" suppresses the singularity. To make this apparent we use the 
completeness of the physical Hilbert space. The physical Hilbert space is spanned by the 
free states. Therefore we can write the identity operator as the sum over all states: 

i=5:/d=^pM°)(p))(«(")(p)i. (72) 

We use this abstract notation because we don't want to bother with conventions which 
are not relevant. The normalization follows from the idempotency of the identity operator 
(1 ■ 1 = 1). 

{u^°'\p)\u'^^\p')) = u^°'\p)u'^^\p'){p\p') = d'^H^p-p'). (73) 



We can project 7+/ (2p+) onto the physical states by applying the identity operator on 
both sides: 



^P al3 •' 

= E/ d'p'd'p\u^''\p'W^^-^^^^{u^^\p)\ = ^ (74) 

There is no mixing between upper and lower components because they are spectrally 
separated. 

Wherever 7^/(2p^) appears we can replace it by l/(2m), since eq. ([7^ ) is an operator 
identity on our space. The instantaneous interaction can be interpreted as nothing but a 
point interaction, 1/ (2p^) being the phase space that goes with it and 7"*' the vertex. 

However, these arguments do not hold in a Feynman diagram. The spin plays a 
dynamical role. In contrast to the instant-form dynamics, where all components of the 
angular momentum are kinematical, in front-form dynamics only the z-component of the 
spin is a kinematical operator. The other components are involved in the interactions. It 
turns out that we can combine the 7"'"/(2p"'") singularity with the l/p"*" singularity that 
appears in a corresponding l.c.t.-ordered diagram with a propagating fermion line, such 
that the singularities of the two cancel. Thus we find that in a Hamiltonian approach 
to light-front field theory, there is an intimate interplay between the singularities of the 
propagators and a singular piece of the interaction. Interestingly enough, this piece occurs 
even for free particles. In old-fashioned ordinary time-ordered perturbation theory, i.e., 
instant form dynamics, the amplitude of propagation depends on the off-shell energy only, 
not on the polarization. 



5.4 Analyticity and covariance 

A Feynman amplitude is an analytical function of scalar objects like pf,Pi- Pj. Often, 
the real values of the scalars are the boundary values of the complex domain on which 
this function is defined. We use these arguments in order to be able to apply the residue 
theorem to integration over the loop variables and to perform Fourier transformations. 
All contour integrals are finite (if we don't pinch the contour) and coincide for integrals 
convergent on the real axis with the integral along this axis. (A coordinate transformation 
y = x~^, as used in the section on zero modes, doesn't alter the results, since it maps the 
real axis on the real axis.) We cannot use the exponential e^^^ to improve the convergent 
along the semicircle of the contour (Jordan's lemma). 

If we integrate over one coordinate separately, manifest covariance is lost. In the case 
of integration over ordinary energy this was not much of a problem since we can consider 
p'^ + m'^ as real and then analyticity in is directly related to analyticity in p^p'^. In light- 
front coordinates the situation is more complicated since p^ — = 2p~p'^ — [p^ + rn^). 
For real values of the scalar object the complex values of p~ and p^ are related {{p\ + m'^) 
is real). The coordinates are each others complex conjugate. After Wick-rotation (p° — * 
—ip^) this remains almost exactly true: p'^ p, p^ ^ —p. For a strip along the real 
axis we have: p~Qp'^ = —p~^Qp~ . Singularities that occur in a complex function can be 



regularized, but the relation between the conjugate variables restricts the possibilities of 
regularization. Singularities of an integer order (like l/x) cannot be integrated by parts. 
But one can approach these singularities in parametric space x"; a —>■ —1. The advantage 
of this dimensional regularization is that it does not interfere with algebraic rules; the 
regularized distributions satisfy the same relations as the singular ones, which is of great 
important for complex, analytical functions. For regularization of complex distributions 
one subtracts these poles as function of the order but with a fixed difference between the 
order of the singularity of p with respect to p: p°'p"-+^ with a fixed k fl^ . 

We will follow a simpler approach with the same result. To avoid complications we 
define distributions of covariant objects only. Analyticity of the covariant object tells us 
the relation between p"*" and p~ at regularization. A homogeneous distribution is given 
by the partial integration of a singular, but integrable, function: 




dp —d) = I dp 

p+ 



ln(p^p ) 



dp 



(76) 



We need p^p^ for positive imaginary values. So we take the cut of the logarithm along 
the negative imaginary axis, therefore the logarithm has an imaginary part of the form 
i7r6{—p'^p~). The homogeneous distribution is 

-(In \p+p-\ + tne{-p+p-)) = — . , _, = PV— - t7Ta{p-)6{p^). (77) 



dp^ p+ + iea{p~) p^ 

The result is that dp+dp- \n{p^p~) = dp- {p^)~^ = 21116 {p^)6{p~^), which is nothing but the 
Mandelstam-Leibrandt regularization. (The i can be accounted for through Wick rotating 
the variable.) After Fourier transformation the distribution \n{p~^p~ + ie) becomes a 
singular function which contains the intersection of the light cone with the null-planes 
x+ = V a;- = 0: T[\n{p+p~ + le)] = = -S\x)5\x+x-) + i6\x±)FV^. This 

fact is part of the reason why there exists confusion about the instantaneous term in the 
fermion propagator. In a Feynman diagram the integrands are treated as meromorphic 
functions in the complex plane. The real part is automatically complemented with a 
imaginary part. If we use a Hamiltonian approach we have constraints which relate real 
parts to real parts and we express one in terms of the other, therefore we might loose some 
information concerning the behavior of the imaginary parts. To put it in other words: the 
off-shell behavior comes naturally in a Feynman diagram and this does not always happen 
in a Let. -ordered diagram. This is another reason why we chose to define light-front field 
theory in close connection with a formulation in terms of Feynman diagrams. 



6 Proof of equivalence 

In this section we have collected several technical aspects of the proof of equivalence. In 
sect. H we gave some examples, in order to illustrate the general procedure. Therefore, 
we do not give any examples here. 

6.1 Energy integration 

We present here the proof of the basic theorem on the integration over p~ . First we 
discuss the case of a single loop. The proof for several loops comes next. 



6.1.1 One loop 



The integration of a Feynman diagram over one energy loop variable p gives the following 
expression the ^'''-interval corresponding to {{QHi > A i < m) V C^Hi < A i > m)}: 



FD{H) 



dp 



27r 2>+ 



2^(p+---p+)A(iJi,---,if^) 



■Pn[v~ - Hi] 



m 
TjN-2 

N-2 
N 



[p- 



H 



N\ 



H 



Hf 


Hi 





1 


HI 


H2 





1 


HI 


H-in 





1 


tt2 


Hm+1 


1 





Hn 


Hn 


1 






(7S 



The last factor is a complicated mixed symmetric polynomial in the ifj's that we denote 
by {—l)"'Wm,n{Hi, ■ ■ ■ , Hm\Hm+i, ■ ■ ■ , Hn+m)- (The phase factor (—1)^+^ is introduced 
to simplify the final expressions.) A is the Vandermonde determinant {n > 2): 



Mxi 



det[Aij-] 

n—j 



(79) 
(80) 



A well known result is 



A(xi 



n— l,n 

n 



See for instance MacDonald or Fulton & Harris |]13| for properties of the Vandermonde 
determinant. 

Proof 

Depending on the (fixed) values of the (kinematical) p^'s some poles are above and 
some are below the real axis. The integral is computed as 27ri times the sum of the 
residues. The residue of pole p~ = Hi can be written as 



(-1) 



i+i '^{Hi, ■ ■ ■ , Hi_i, Hi^i, 



Hn) 



2^ipt---p%)AiHir--,HN) 
We can add a line and a column to the determinant in the numerator: 



(82) 
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(83) 

The final formula is obtained by adding determinants of type eq. ( p3D which amounts to 
just adding their last columns. This gives the result stated in the theorem. 



6.1.2 Several loops 



In the case of several loops we integrate loop by loop. We must use the residue theorem 
such that the order of integration does not change the result. In general, the momenta of 
the particles on the internal lines will be linear combinations of the integration variables, 
say ki = Y^k<^iPk, where kf is the four momentum on line i and is the integration 
variable. One has the freedom to choose the latter such that ctf is either +1, —1 or 0. 
Theorem Multi-dimensional energy integration 
An unambiguous expression is: 



/n 
dprdP2 • • • dp- n(E <^'p-k+Hi)-' = 
i=l 



\JliJ2i iJmS 



IJl-J 




The antisymmetrized product [a^ ■ ■ ■ Oi"^]ji-j^ is the determinant of the matrix 



i^)jl—jm — ^ 



/< ••• a] \ 



(85) 



The inverse of the matrix a will be denoted by a. The poles {ji, ■ ■ ■ , jm} that are included 
in the sum have for all values of the correct imaginary sign of H.j/a'^.j. because this sign 
is determined by the plus-components of the integration variables. Before we start to 
integrate we first determine inside which contours the different poles lie, then drop the 
ie-description. This is clearly an invariant formulation, so the order of integration can be 
altered. 
Proof 

An unambiguous way to define the integration is to shift the integration contour 
slightly into the complex plane and leave the poles on the real axis. The poles are 
determined by the following set of linear equations 

J2alP7 + Hj, = 0, t,k = l,---m, (86) 

i 

which have the solution 

PT-~j:^'Hk. (87) 

k 

Next the multidimensional integral is written in terms of these new variables. The Ja- 
cobian of the transformation is det« = 1/ det« = ■ ■ ' (^"^]ji---jm- Subsequently, we 

apply the residue theorem to every p~-integration. We substitute the new variables, and 
find that the pole part of the integral factorizes: 

f " f{z) - f 1 f 1 

/ n ^^in ^ — n ^ ^ ^^^E"'] / — ~j — / — ~j — ^^""y"!- 

J j^i {aizj - 'yi) ■ ■ ■ {ainZj - 'Jm) Z/i-ai7j •' Um-am'yj 

(88) 



We used the summation convention. Note that the integral is independent of the choice 
of integration variables. 

This type of multi-dimensional complex integration is not related to topology, so the 
deformation of the contour might change the result. The torus obtained by closing the 
different contours depends on the choice of coordinates. To avoid ambiguities we take an 
algebraic view instead. 

6.2 Recursion formula 

The recursion formula is the basis of the proof of equivalence. It tells us how to take 
out of any Feynman diagram the building block of a l.c.t-ordered diagram: an energy 
denominator. This happens without changing the structure of the algebraic form of 
the reduced Feynman diagram, so we can apply this formula, a number of times (the 
recursion). The final result, obtained upon the last application of the recursion formula 
can immediately be evaluated. The final object is just a piece of a Let. -ordered diagram 
(TOD): a product of energy denominators. 

The recursion formula allows us to consecutively pull energy denominators out of 
FD{H) in order to obtain a sum of TOD^s. 

Theorem The following identity is true for any m and n {N = m + n): 

A{Hi,---,Hn) Hi — Hm+i 

A{H2,---,Hn) A{Hi, ■ ■ ■ , Hm, Hm+2, ■ ■ ■ , H]^) J 

Remark 

The reduction step removes two poles. Hi and Hm+i, and combines them into a single 
energy denominator. Hi — H^+i- The second factor of the r.h.s. of eq. (^) consists 
of two terms, both of which contain one pole less than the original form. The factor 
{Hi —Hm+i)^^ is incorporated in the TOD. At the first stage, Wm,n/ ^ is a structure that 
is directly related to a Feynman diagram. After taking steps in the reduction algorithm 
objects with the same structure are obtained, but these objects are not in the same way 
associated with (possibly different) Feynman diagrams. The last step in the algorithm is 
given by: 

Wi^n{y\Xl, ■■■,Xn) = Wn,l{Xi, ■ ■ " ,X„|y) = A(xi, ■ ■ ■ , X„) (90) 

First we prove the formula (^9]) and then we show how to carry out the reduction. 
Proof 

First express Wm,n in terms of a determinant as in eq. ([781). Then perform the usual 
manipulations with determinants: take linear combinations of rows or columns. If we 
subtract the first row from the rows 2 to m, the m -|- 1'^* row from the other rows {m -|- 
2, ■ ■ ■ ,m + n}, and expand the determinant with respect to the last two columns, we 



obtain: 
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where K — n + m 
obtain: 



2 = 



• yn-yi 

2. We can add a row and a column to the determinant to 
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• yn -yi 

We spht the determinant into two parts by adding and subtracting in the last column the 
column {x2 — Xi, ■ ■ ■ , — 0, ■ • • , 0). Then we subtract the last column from the next 
to last column in both determinants to obtain 
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We add the m-th row to the rows above it in the first determinant, and subtract it from 
the rows below it in the second determinant. The result is 
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= n + m — 3 and define for any k the symmetric function by the relation 
= {x — y)(f)''~^{x, y). The rows contain one of the factors Xi — Xi, yj — Xi, yj — yi 
Xl, that can be divided out. The product of these factors are written as the ratio 



of two Vandermonde determinants, to obtain from eq. (94): 



xi-yi 



A(a;i ■■Xm,y2 ■■yn,yi) 
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A(X2 ■■yn,Xl) 
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(f)'^{Xm,Xi) 

^^Hvuxi) 

(t>^^{y2,Xi) 



1 

1 

1 
1 



(t^'^'iyn^Xi) "10/ 
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The dependence of the first determinant on yi is only apparent. The same is true for 
the dependence on xi of the second determinant. We can easily see this through matrix 
multiplication: 
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V </>"-2(x„,y) 



xi+y 1 1\ 

Xm+y 1 1 
XjnH+y 1 



Xn+y 1 0/ 
(96) 

The determinant of the second matrix at the l.h.s. is 1, so the determinant of the first 
matrix at the l.h.s. is the same as the determinant of the matrix at the r.h.s.. Removing 
the Ui dependence in the first determinant and the Xi dependence in the second we get 
two familiar objects: 
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These determinants are nothing but ly-functions, but now with less arguments than we 
started with. So we can write eq. ( p7D as follows 

A(a;i ■■■yn) / Wm ,^-1(3^1 ■ ■ ■ Xm\y2 ■ ■ ■ Z/n) Wm-l,n{x2 ' ' ' ^mlVl ' ' ' 2/n) | (98) 

If we divide the whole expression by A{xi ■ ■ ■yn) we get the reduction formula. 



6.3 Reduction algorithm 

We will now show that the application of the formula derived above gives us parts of 
l.c.t-ordered diagrams. First we need to specify the structure of a Let. -ordered diagram. 
Definition Loop connection tuple 

The loop connection tuple H = {Hi, ■ ■ ■ , iJ„) is an ordered set of objects related to 

the propagator denominators — + ze = 2pj{p' — Hj) or Hj = p~ — pj + + ^ ■ 



The ordering of the tuple corresponds to consecutive ordering of the internal lines in the 
corresponding loop in a Feynman diagram. 

We will use the terminology of "lines" when we mean the corresponding momentum 
or energy, or state. There is some arbitrariness in the definition of the momenta in the 
loop, but in the objects of interest Hi — Hj this arbitrariness is gone because they are 
invariant under a shift of the loop momentum (p^ — p'^ + a^) . The expression — Hj is 
the total incoming J2 Pext momentum minus the on-shell values of the minus-momenta of 
the internal lines P~ on- shell + Plan- sheib calculated with the help of p+- and pf^^. (See also 
sect. ^) 

Definition Backward and forward 

A line i of the loop connection tuple is going backward if the object Hi has a positive 
imaginary part and is going forward if it has negative imaginary part. Thus in the 
Feynman diagram above, eq. (^), Hi, - ■ ■ , Hm are going forward and Hm+i, ■ ■ ■ , Hjsi are 
going backward. The sign of the imaginary part is opposite to the sign of the on-shell 
energy of the particle, therefore this definition of backward and forward coincides with the 
causality condition: positive-energy particles go forward in I.e. time and negative-energy 
particles go backward in I.e. time. 
Definition Early, late and trivial events 

An early event is a vertex between a backward and a forward going line, a late event is 
a vertex between a forward and a backward going line, if one goes around the loop in the 
order corresponding to the connection tuple. All other vertices are trivial events. There 
are equal numbers of early and late events. 
Definition Flat loop diagrams 

A flat loop diagram has one early (and thus one late) event. 
Definition Crossed loop diagrams 

A crossed loop has more than one early event. 
Definition Skeleton graph 

A skeleton graph is the tuple of signs of the objects Hi of a connection tuple. It is 
given by the mapping {Hi, H^, ■ ■ ■ , Hn} — > {crC^Hi), a{QH2), ■ ■ , a{QHn)}- The function 
a is the sign function. 

For different external momenta in the Feynman graph we have a different set of skeleton 
graphs, and for each Feynman diagram there are a number of skeleton graphs. For each 
sector (associated with a specific number of poles above and below the real axis) of the 
loop momentum p'^ there is a skeleton graph, thus for a loop with n lines there are n — 1 
skeleton graphs. 

The skeleton graph already tells us the general features of the l.c.t.-ordered diagrams 
which are contained in a Feynman diagram, because it tells us the direction of the internal 
lines. This is used as our guide how to take "time-slices" of the Feynman diagram. The 
direction of a line tells us in which order we can encounter events (vertices). Early and 
late events correspond to sign changes in the skeleton graph. 
Definition Causally connected events 

Two events are causally connected if they lie on the same loop and there are neither 
early nor late events lying in between. 

So, two causally connected events are connected by parts of a loop that consist of lines 
that are either all forward or all backward. 
Remark 

Clearly, it makes sense to say that causally connected events are ordered in I.e. time. 



If we follow a loop in the direction given by the orientation defined by its connection tuple, 
then we will encounter forward and backward going lines. If two vertices are connected 
by a forward line, they are said to be ordered in I.e. time in the same way as they are 
ordered in the loop. Otherwise their order in I.e. time is the inverse of their order in 
the loop. This partial ordering, which is given by the skeleton graph, is obviously not 
complete. Only causally connected events are mutually ordered this way, but not with 
respect to other events. 

Note that we don't make any statements here about reducible Feynman diagrams, 
which is a completely different story. Our causally unconnected parts connect up a later 
time, so they are parts of the same irreducible Feynman diagram. 
Definition Simultaneous 

Two parts of a skeleton graph are said to be simultaneous if they do not share events 
that are causally connected. 
Remark 

The flat box that we discussed in sect. ^ consisted of a late and an early vertex, 
connected by two distinct parts of the loop, one consisting of lines graded +, the other 
of lines graded — . The relative l.c.t.-ordering of the events on these two parts is not 
necessarily determined by the skeleton graph, but application of the reduction formula, 
eq. (^91), produced immediately the two possible Let. orderings. The diagrams found 
showed the expected energy denominators. 

In situations where two simultaneous parts occur, the reduction formula does not 
provide immediately the l.c.t.-ordered diagrams. An example was given in sect. |2.1.2| . 
From the point of view of I.e. time ordering, one expects diagrams to occur corresponding 
to all relative Let. orderings of simultaneous parts. In momentum-energy language this 
means that diagrams with certain energy denominators should occur. Indeed, this is the 
content of the next theorem. 

Theorem Simultaneous parts come in all combinations. 

Remark 

The proof of this theorem relies again on a recursion. First we suppose that we have 
two simultaneous parts, that are already ordered by themselves, but not mutually. Both 
parts correspond to sets of energy denominators, say {a} and {/?}. So we have the two 
distinct TOD^s Ila^"^ and 11/5^"^. In this language the content of the theorem can be 
written as follows, {{ik, jk}k\[{ik, jk} ^ {ik+i,jk+i}] A [k < I ^ (ik < k) A {jk < ji)]] 

n«r^n/?7'= E nK+/?.j-'- (99) 

i=l j=l all ik,jk k=l 

Proof 

We apply the formula: 



i=i j=i 



X 



m— 1 n m n—1 

i=l 7 = 1 i=l 7 = 1 



(100) 



We can apply this algorithm recursively to obtain TOD^s with energy denominators 

of the form Q^i + S Pj- If the product consists of more than two TOD^s we can apply 
this algorithm recursively, to two TOD^s at a time, since it is associative. 
The expressions obtained are composite energy denominators: 



(p-(a) - Ho{a)) + {p-{(3) - HoiP)) = P-{a U (3) - Ho{a U (3). (101) 



6.4 Reduction of Feynman diagrams 



The previous sections were mainly concerned with algebraic identities. Now we turn to the 
general strategy of the reduction. The reduction algorithm for a flat loop is straightforward 
as we noted when we discussed simultaneous parts. The trivial vertices come in all 
orderings of vertices on the forward line C^H < 0) with respect to those on the backward 
line {QH > 0). The vertices on one hne are already ordered with respect to each other 
by the skeleton graph. Starting with the early event one can reduce the lines next to the 
early event. This algorithm ends and gives 2^~^ TOD^s if the loop is flat for all skeleton 
graphs. For each skeleton graph there are (m-i) TOD^s. 



The mapping from a recursive algorithm to a time-ordered theory is straightforward 
for the fiat loop. The order in which the poles Hi are removed in the reduction formula 
is the same as the time ordering. We refer to this reduction as time ordered reduction. 

A flat Feynman diagram gives the topology of all TOD^s contained in it and all com- 
binations between trivial events appear respecting the causal order. Conservation of p'^- 
momentum determines whether a trivial vertex is an absorptive or emissive event. How- 
ever, this is not a new element, because we have seen that a skeleton graph is determined 
by the external lines besides the value of the plus-component of the loop momentum. Of 
course, early events absorb external particles while late events emit them. 

Because Feynman diagrams form the basis of our treatment of light-front field theory, 
the basic elements we are concerned with are the single particle propagators and the 
vertices derived from the underlying Lagrangian. The interaction Hint is derived from Ci^t 



with an additional phase factor yy'^'^Pipt • • -Pn j where the longitudinal momenta 

of all the lines from a vertex are included. A wave function must also be multiplied with 
the phase factor, as compared to the covariant wave function. 

For the success of our reduction algorithm the details of Cint are of minor importance. 
The algebra is connected to plus-momentum flow in loops. The only internal lines in the 
diagram of interest are those in the loop, how momentum is extracted from and added to 
the loop is of less importance. 

The most important feature of the reduction algorithm is the fact that it always starts 
from an early state with positive longitudinal momentum (there is always an early state 
as the result of momentum conservation), so we exclude "vacuum" -type diagrams. 

First, reduction is performed on the skeleton graph starting from the early events, i.e., 
removing lines directly connected to the early events. This can be followed by removing 
poles corresponding to consecutive pieces of the loop until a late vertex is reached. This 
is the point where Wm,n/^ is reduced to a sum of terms of the form Wj^i/A' or W^i^^/A". 
Secondly, the simultaneous-parts theorem is applied to write all these terms as sums of 
terms containing true energy denominators 

Because we could have started the reduction from the late vertices in stead of the 
early ones, we see that the algorithm can be written in different ways. The final result, 
however, is the same. The same is true, mutatis mutandis, for the application of the 
simultaneous-parts theorem. 

So we peel off more and more of the Feynman diagram in a manner which is locally 
(for causally connected events) equal to time ordered reduction. 

The relation of the results of the reduction process to l.c.t.-ordered diagrams in the 
crossed loop case is more complicated than in the case of a flat loop. The simple heuristic 



Mapping 




of cutting lines representing constant I.e. time surfaces leads to a more elaborate book- 
keeping in the crossed-loop case. This is so, because only the global structure of the 
Feynman diagram determines which simultaneous parts are joined by early or late ver- 
tices. The general strategy proposed here is to start at an early vertex, use the reduction 
algorithm locally until a late vertex is attained. This procedure is to be repeated until all 
late vertices have been processed. Next apply the simultaneous-parts theorem repeatedly. 
Extension mapping 

Multi-loop diagrams can be reduced one loop after another. The loop momenta that 
are not integrated over are kept fixed. The skeleton graph tells us again what is the 
general form of the time-ordered diagrams. Therefore, the skeleton graphs, associated 
with different finite domains in p+-space need to be determined first. 

Upon application of the reduction algorithm to the first loop, energy denominators 
occur that are combinations of two propagator poles. When the next loop is treated, 
some poles come from those energy denominators, while the others are due to propagator 
poles occurring in the part of the original Feynman diagram unaltered by the reduction 
so far. These different types of poles play the same role in the integration of the next 
variable. The pole is again the difference of the p~-flow and the on-shell energies from 
the poles and their associated lines in the Feynman diagram. The question which lines 
must be combined to generate energy denominators is related to the imaginary parts of 
the propagator poles and thus answered when the skeleton graph is determined. During 
the reduction process these combinations remain fixed. 

The integration is invariant under coordinate transformations and reordering of the 
integrations. The mapping from a recursive algorithm to a l.c.t.-ordered approach is 
more complicated here than in the single-loop case, but in essence the same. The most 
complicated task is the construction of the skeleton graphs. After this job is done, the 
reduction algorithm is applied to one loop after another, and the interpretation of the 
result is the same as in the case of one loop. 
Theorem Spectrum condition 

The spectrum condition p+ > holds for all particles in the internal loops of the 
Feynman diagram. 
Proof 

The spectrum condition follows from two ingredients. First, at any vertex there is con- 
servation of four-momentum, in particular plus-momentum. Secondly, lines with negative 
p"*", antiparticle lines, can be reinterpreted as particle lines with positive p+ by reversing 
the direction of the four momentum on such lines. This reversal is in agreement with 
four-momentum conservation and Let. ordering. 



6.5 Algorithm proofs 

In a number of cases we restricted ourselves to a recursive formula, which could be applied 



to any structure [Tj]. We gave a rule that could be applied successively and leads to l.c.t- 
ordered diagrams with the right structure. For the flat loop case we obtain all l.c.t.-ordered 
diagrams, since the algorithm is unique, and directly related to a time-ordered picture. In 
the case of the crossed loop diagrams the algorithm is not unique; the rule could be applied 
to different parts of the structure at the same time. We could start with the reduction 
at one early vertex or another early vertex. Strictly speaking, this poses a problem only 
when the rule can be applied to overlapping regions of the structure, the different orders in 



which the rule is apphed might lead to different answers. Our conjecture is that in our case 
the final answer is unique and does not depend on the order in which the rule is applied. 
We base this conjecture on the observation that each term obtained after applying the 
first rule (with the simultaneous parts still present), can be associated with a topological 
object, such that all terms give the complete set of objects with specific properties. The 
topology is invariant, thus unique and independent of the way it is obtained. We have 
checked this for a number of cases, and found it to hold in all these cases. Note that we 
used a similar argument for deriving an invariant multi-loop integration. 

We would rather have shown that the reduction algorithm would lead to all time- 
ordered diagrams which satisfy the spectrum condition, which would also be some kind 
of invariant, but this turned out to be too complicated a problem in the general case. 
Therefore we did not write down a general formula, directly relating any Feynman diagram 
to the sum of Let. -ordered diagrams. 

At the moment we satisfy at each step the spectrum condition and allow all possible 
orderings. The global properties follows from the local ones. 



7 Discussion 

We have established the degree of equivalence between light-cone-time ordered pertur- 
bation theory and covariant perturbation theory for spin-0 and spin-1/2 particles. This 
effort might seem superfluous since the connection between ordinary time-ordered pertur- 
bation theory and covariant perturbation theory is well established [T^. One might be 
tempted to believe that the methods that work in the case of ordinary time- dependent 
perturbation theory apply to the Let. -ordered theory too. This belief belongs to folklore. 
In practice the understanding of light-front field theory is growing only slowly. Basic 
results in covariant field theory were often proven to obtain also in light-front field theory 
along a path through much trial and error. It took years to obtain a proper light-front 



version of the Schwinger-model [0 and to prove spontaneous symmetry breaking in the 
light-front version of 0^ fl^. The renormalization of light-front versions of known covari- 



ant, renormalizable field theories is still an unsolved problem |T8[. In general, approaches 
are followed that are specially tuned to the problems of light-front field theory. Therefore 
it turns out to be difficult to relate the solutions obtained to basic results in covariant 
field theory [||. 

Conceptual problems were already present from the onset of light-front field theory, 
when Weinberg showed ||^ that only some physical processes, each represented by an 
ordinary time-ordered diagram, contribute to the Feynman diagram if this diagram was 
calculated in a frame that moves with the speed of light This so called infinite 



momentum frame (IMF) cannot be connected to any other reference reference frame by 
a finite Lorentz transformation. Thus, a limiting procedure is involved. This limit has to 
compete with other limits present in field theory: infinite space integration, regularization 
of singular expressions. This is separate from the additional problems that might arise 
when we are dealing with fermions That the IMF is naturally described with light- 
front coordinates is only apparent for coordinates and momenta That only some 



diagrams survive 0] is puzzling. However, we have shown that this, in general, is indeed 
true. Another possible approach to light-front field theory is the direct quantization on 
a light-front. There are many different ways to do this, which become more elaborate 



if the theories are supposed to incorporate more features [24]. As a classical theory, 
light-front field theory is ill-defined; the standard initial values problem on the light- 
front is overdetermined. In addition, it leads to a non-unique evolution in time ^j. The 
first problem can be solved in principle through methods devised for the quantization of 



constrained systems |2^. The second problem is more serious. One needs to introduce 
degrees of freedoms associated with different evolutions in I.e. time and then introduce 
constraints which can restrict the space of solutions to the one considered physical on 
some grounds. This, in essence, is what people are dealing with when they introduce 
zero-modes, degrees of freedom of which the evolution is unknown (zero or infinite?). 
This problem is often disguised in practical calculations where a ^-singularity occurs 

Another way to quantize a light-front theory is to use axiomatic commutation relations 
on a complete set of free fields an approach guided by the results of current 

algebra |^. On the light front, different points with (Aa;_L = 0) are light-like separated. 
The question arises what is the equal-time commutator between fields: a delta function 
in x~ which violates covariance, or a sign function in x~ which leads to non-integrable 
fields (except at the loss of covariance) |3l|. Whether such theories can describe physical 
processes has not been established. This approach became less favorable in the late 
seventies. The approach most favored nowadays is based on two methods: with the 



covariant results in mind derive a "constrained Hamiltonian" 0, |3^, [^. Due to zero- 
modes it is hard to make a one-to-one correspondence between normalized states on a 
space-like surface and a light-like surface. As long as one deals only with tree graphs 
all these problems are rather formal. The presence of loops makes them acute. In loops 
one has to "sum over all states". This rule forces one to think over states with — > 0. 
The advantages of light-front field theory are paid for by the occurrence of problems. In 
the approaches discussed here the problems arise in different forms and are dealt with in 
different ways. They might be disguised as technical difficulties. Due to the fundamental 
nature of these problems the final results depend strongly on the choices one has to make 
when defining the (finite) theory, e.g., boundary conditions of fields quantized in a box 
3^ , []35| , or regularization of the (p"*")"^ singularity ||2^. We also had to make some 



choices. Wherever we had to do so, we emphasized the relation with Feynman diagrams. 
In a manifestly covariant approach there is no (p"*")"^ singularity. It is a distinct advantage 
of our approach that this singularity also does not occur (see sect. We seem to have 
cured one of the diseases of light-front field theory. However, presently we do not know 
whether our regularization procedure leads to the same answer as the covariant approach. 

We have shown (see sect. pTT| ) that in some cases there are S{p'^) contributions to the 
S-matrix. Maybe these terms indicate a coupling to the "vacuum" or they may represent 
contributions which relate one version of light-front field theory to another by a finite 
renormalization. However, it would be good practice to try to separate the mathematical 
question "how to calculate?", from the physical one "how to interpret?". 

We consider the present situation in light-front field theory to be confused. We give 
three reasons for this point of view: 

(i) The paucity of comparisons to standard covariant theories; 

(ii) The mixing of mathematical problems with physical ones; 



(iii) The lack of consensus on what are the established results (with proper, mathemat- 



ically rigorous derivations). 
Still, there are a four good reasons to work on light-front field theory: 

(a) It is the only theory distinctly different from covariant field theory which allows for 
a comparison at intermediate levels. Such a comparison increases the understanding 
in both theories; 

(b) It is the most natural way to describe nucleons in terms of quarks; 

(c) Our understanding increases with each answer to questions that light-front field 
theory raises; 

(d) It is a Hamiltonian formulation, sharing the intuitive picture inherent in elementary 
quantum mechanics. 
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